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Heavy hadrons on Nf = 2 and 2 + 1 improved clover-Wilson lattices 

Tommy Burch^ * 

^D-93053 Regensburg, Germany 

We present the masses of singly {B, Bs, A;,, Ej,, etc.), donbly {Be, gs, T, Hi,c, 'S.bb, etc.), and 
triply {Q.bcc, ^bbe, i^bbb, ctc.) heavy hadrons arising from (QCDSF-UKQCD) lattices with improved 
clover-Wilson light qnarks. For the bottom quark, we use an 0(a, i!"^)-improved version of lattice 
NRQCD. Part of the bottomonia spectrum is used to provide an alternative scale and to determine 
the physical qnark mass and radiative corrections used in the heavy-qnark action. Results for spin 
splittings, opposite parities, and, in some cases, excited states are presented. Higher lying states 
and baryons with two light quarks appear to be especially affected by the relatively small volumes 
of this (initially) initial study. This and other systematics are briefly discussed. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.39.Hg, 14.20.Mr, 14.40.Nd 


I. INTRODUCTION 

Ditto the above [T]. 

Here are some references concerning experimental and 
lattice-QCD results for heavy-hadron spectroscopy: 

New bottomonia, including the ground-state rjb [SHI] 
and the corresponding “radial” excitation r]b{2S) |Sl[n|, 
the spin-singlet P-waves hi,{lP) and hb{2P) [7], and a 
possible Xbi^P) |S|; 

Orbitally excited B and Bg mesons (SHU); 

4*^ na, Sf, [Him, and Qb baryons [TMIB] - Excited 
Sfc m and Ab baryons [snuni; 

Lattice studies of bottomonia , including those 

with charm sea-quarks [53 123: predictions of D- 

wave [27] and higher states [28] - 

Lattice results for B and Bg mesons [29H36] , including 
predictions for the B^ system mm-. 

Lattice results for 6-baryons [3ll |39l [40] and triply bot¬ 
tom baryons mm- 

Some more recent developments: Lattice studies of be 
baryons [33] i bottomonia [33] i positive parity Bg 
mesons [45| : and the experimentalists have been busy 
with the Xb{3P) state and Sf, baryons [36H38] . 

For heavy-quark-model calculations, see, e.g.. Ref. [39] . 

Once again [T]. 


II. LATTICE CALCULATION 

In the present section we give the details of the sim¬ 
ulations: the gauge configurations used, how the quark 
propagators are calculated, and how these are put to¬ 
gether to form the correlators for the hadrons of interest. 


A. Configurations 

We use gauge configurations which include either Nf = 
2 flavors of non-perturbatively improved clover-Wilson 
quarks [^ or = 2 + 1 flavors of SLiNC quarks m- 
The relevant parameters of the ensembles used can be 
found in Table |l] As can be seen in the last column, the 
spatial extent of the lattices is rather small. Whereas this 
may not strongly affect tightly bound, multiply heavy 
systems (e.g., T, Be, i^bcc, etc.), higher excitations and 
hadrons with lighter valence quarks may “feel the pinch” 
[55] and this is therefore a source of systematic error 
which we must keep in mind. 


B. Light and charm quark propagators 

The propagators for the light {u,d), strange (s), and 
charm (c) quarks were produced using the Chroma soft¬ 
ware library [53] . The light-quark propagators on ensem¬ 
ble a and the light- and strange-quark propagators on en¬ 
sembles b and e were originally created for other projects 
(see [50] and [53] i respectively) and we must work with 
the (rather severe) quark-source smearings chosen therein 
(see Table |nj of course, with great computational ad¬ 
vantage of just having to read in the files). The light- 
and strange-quark propagators on ensembles c and d 
and the charm-quark propagators on ensembles b-e were 
created with the aim of better resolving excited states 
(less smearing, see Table [h]). The charm-quark mass was 


TABLE I: Relevant lattice parameters [501151] . 


Ibl 

/5 

N'^ X Nt 

TV 

aM,, 

M^L 

a 

5.29 0.13632, - 

32^ X 64 

660 

0.1070(5) 

3.42 

b 

5.50 0.12104, 0.12062 

24^ X 48 

320 

0.1406(8) 

3.37 

c 

5.50 0.12100, 0.12070 

24^ X 48 

180 

0.1515(10) 

3.64 

d 

5.50 0.12095, 0.12080 

24^ X 48 

210 

0.1661(8) 

3.99 

e 

5.50 0.12090, 0.12090 

24^ X 48 

800 

0.1779(6) 

4.27 
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taken from a related study [Ml (kc = 0.1109). The 
lattice scale there, however, was set using a different ob¬ 
servable (the flavor-singlet baryon-mass combination Xjq 
[5T] ) than the one here {M{1P) — M(1S) from bh) and we 
therefore have a systematic shift in our i?c-system masses 
when using the bottomonia scale (see Sec. IIIB). 


C. Non-relativistic quark propagators 


For the bottom quark, we employ improved NRQCD 
[59]) including terms up to 0 (u‘^), where v is the heavy- 
quark velocity. We use the time-step symmetric form of 
the evolution equation: 

ai/o(t + a)\" 

2^1 ) 

y _ 

•(()(x,t), (1) 


(j){y,t + a) = 


^ _ aSH{t + ay ^ 




where the binomial expression of the exponential of the 
lower-order (in v) terms is carried out to n = 4. Hq han¬ 
dles the heavy-quark kinetic energy and the associated 
0(a) time-step correction. 


Hn = --. 


A 


A2 


2mQ An 4toq ’ 
and 6H contains the O(u^) relativistic corrections. 


SH = -^A2 


8m| 

igc^ 

8ml 

gcs 

8 w?q 

gcA 

2mQ 


(V-E-E-V) 
cr • (V X E - E X V) 


cr-B . 


( 2 ) 


(3) 


TABLE 11: Quark source smearings. 


Ibl quark 

sm.type 

params. 

a u,d 

Gauss l^l -1- 

(k = 0.25, A = 400)+ 


APE [58] 

(/ = 2.5,A = 25) 

b u,d,s 

Gauss -1- 

{k = 0.25, A = 150)+ 


APE 

(/ = 2,A = 20) 

c,d u, d, s 

Gauss -1- 

(k = 0.25, A = 20)+ 


APE 

(/ = 2,A = 3) 

e u,d,s 

Gauss -1- 

(k = 0.3, A = 130)+ 


APE 

(/ = 2,A = 20) 

b-e c 

Gauss -1- 

(k = 0.25, A = 12)+ 


APE 

(/ = 2,A = 3) 

a-e Q 

Gauss 

(k = 0.2, A = 16) 


Covariant derivatives and Laplacians with a tilde rep¬ 
resent improved versions, where next-to-nearest neigh¬ 
bor sites are included (this is the same form of correc¬ 
tions we used in previous studies of bottomonia ED); the 
form of the heavy-quark evolution, Eq. Q, has been im¬ 
proved for the present study m)- The (chromo) elec¬ 
tric and magnetic fields are constructed using the stan¬ 
dard four-plaquette clover formalism (see, e.g.. Ref. [52]) 
and tadpole improvement [MIES] is applied throughout: 
Uf^{x) Uf_,{x)/uo, E(a;) E(a;)/M^, B(a:) B{x)/uy 

where uq is determined from the average plaquette. On 
each ensemble we run two or three different heavy-quark 
masses {amg = 1.5 — 3.0) and for the radiative correc¬ 
tions, we use tree-level values Ci = 1 , as well as the case 
where C 4 = 1 . 2 . 


D. Hadron correlators 

As already mentioned above (Table[n| , we use a combi¬ 
nation of source smearings for all quark propagators. For 
the heavy quark (Q), we use the smeared source, as well 
as another where a covariant Laplacian is also applied 
(giving a radial node; for P- and H-wave mesons, we use 
a local source as the second choice). For all quarks we use 
local sinks, while for the heavy quark we also consider 
the two smearings used at the source. This leads to a 
rather peculiar situation where the heavy-quarkonia and 
triply-heavy-baryon correlators form 2x3 matrices, with 
a 2 X 2 symmetric block, whereas all correlators involv¬ 
ing light, strange, or charm quarks, together with heavy 
ones, form 2x3 off-diagonal blocks of a larger (mostly 
unknown) matrix. This is a not a major problem, how¬ 
ever, as we can still fit such heavy-light correlators to the 
usual ansatz, 

C(% = (0|0.(t)0](0)|0) 

00 

= (4) 

n—1 

except that we must fix one amplitude for each energy 
level considered (the fit then gives amplitude ratios, but 
the same energies). For some fits, we find it advanta¬ 
geous to consider a submatrix (2 x 2 or 2 x 1 ) of the 
ones we have (this is likely due to limited statistics) and 
for most mass differences reported, we use appropriate 
(jackknifed) combinations of only the smeared-source, 
smeared-sink correlators (see below). 

The interpolating operators that we use to combine the 
quarks together into the mesons of interest are shown in 
Tables HTTl and Hvl One needs to be careful when combin¬ 
ing the u, d, s, c quark propagators from Chroma with the 
nonrelativistic Q propagators: a change in the spin-basis 
is needed [M| . 

The baryon operators can be found in Table|V| In order 
to project out the desired spin and parity, the baryon 
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correlators should then be of the form 

S'/"^(i) = (^(l±74)00^ (5) 

or, for operators with an open Lorentz index, 

= , ( 6 ) 

where the zero-momentum spin-projectors are = 

Sik - \lilk and = | 7 i 7 /c. In the end, we average 
over the nine remaining spatial indices For the 

flavor projections, we follow lowest-order HQET and do 
not consider mixings between the different heavy-quark 
configurations (e.g., between Aq and Ag in Table jv] we 
mostly consider the heavy-light diquark configurations in 
order to reach the negative-parity states). Although it is 
a poorer approximation, we do the same when consid¬ 
ering different charm-quark configurations {q or q' = c). 
Depending on the desired state, however, it may be that 
one must take care with the appropriate flavor projec¬ 
tions of the light quarks (g, q' = u, d, s) [S5] . 

We also create non-zero-momentum correlators 
(smeared source / local sink only; p = 2'Kn/L, where 


TABLE III: Quarkonia operators. 


lowest state 

irrep A 

operator 

Vb 

0 -+ 

Ai 

xV 

T 

1 ““ 

Ti 

x'^i<l> 

XbO 

0 ++ 

Ai 


Xbl 

1 ++ 

Ti 

X} kfj) 

Xb2 

2 ++ 

T2,E 

CfcVfe)!/) 

hb 

1 +- 

Ti 

x^Vib 

Vb2 

2 "+ 

E 

^ j'^ j)4’ 

T2 

2 ““ 

E 

X^(ViVjCrfc -I- VjVk(yi)4> 


TABLE IV: Heavy-light meson operators (for the relativistic 
quark: q — (“) [M|L 


lowest state 

^ min 

irrep A operator 


0 " 

Ai 

yW 

b; 

1 " 

Ti 

X^(JiU 

^qO 

0 + 

Ai 

X^l 

^ql 

1 + 

Ti 

X^CFil 


0.2 — 


CL, 


0.1 — 


32®,N,=2 + 
am(,= 1.5: 24®,Nf=2+l □ ; 32®,Nf=2 K 


0.0 


0.00 


0.01 


0.02 
(a MJ® 


0.03 


0.04 


FIG. 1: Spin-averaged IP — IS QQ mass differences versus 
pion mass. 


|n| < 3) for the T, B*, and P* in order to determine 
the kinetic masses of these mesons from their dispersion 
relations. This provides us with an absolute mass scale 
and a way to set (or interpolate to) the physical 6-quark 
mass. 

For many of the heavy baryons considered herein, we 
present an alternative analysis in which we try remove 
most of the remaining, leading uncertainties in the heavy- 
quark parameters. By considering appropriate combi¬ 
nations of (jackknifed) average correlators, we subtract 
E{T)/2 for each 6 quark and E{Bc) — E{T)/2 for each 
c quark. Inserting experimental values for M(T) and 
M{Bc) into these mass differences, we arrive at more 
precise, absolute estimates of the heavy baryon masses, 
albeit via “less predictive” means. 


III. ANALYSIS 

In the following subsections we present our analysis of 
the heavy hadron correlators, leading to our results for 
the associated masses [55] . 


TABLE V: Heavy baryon operators (for the NR quark: Q = 

(o") m)- 


label 

states 

. 

^ mm 

Sdq 

operator 

Aq 

Ai), —be? ^bc 

1 + 

2 

0 

{qCxbq')Q 

A^q 

a(*) 

^-^bc ^^‘‘bc 

1± 

2 

0 

{qCxiQ)q' 

'Eqi 

yiC*) >*) >*) 

^b ’ “b ’ “b ’ “be ’ “bcc’ •• 

1+ 3 + 
■2 ’2 

1 

{qCxiq')Q 

^Qi 

„(*/*) q(*/*) 

Zjb ? “b ? i^b ’ ••• 

1± 3^ 

2 ’ 2 

1 

{qCxiQ)q' 

SqQz 

■^(*//*) 

“bb ’^^bb ’ “bbc 

1± 3± 
2 ’ 2 

1 

{QCxiQ)q 

^QQQi 

^bbb 

3 + 

2 

1 

{QCxiQ)Q 


A. Quarkonia 

At the present stage in this project we have not yet 
achieved a high-precision analysis of the bottomonia sys¬ 
tem; more precise studies may be found in Refs. [HHIHI- 
However, we need to start somewhere and, as such, we 
are interested in using part of the 66 spectrum to set 
the lattice scale and the parameters used in the heavy- 
quark action (mg, c^). We present a few other results 
(spin splittings and excitations) along the way. Further 
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quarkonia simulations (e.g., on higher-statistics, larger- 
volume Nf = 2-1-1 ensembles), hopefully leading to a 
more complete analysis, are currently underway [T]. 

For spin-averaged and spin-dependent splittings we use 
single-elimination jackknife to create appropriate combi¬ 
nations (e.g., ratios or ratios of products) of the smeared- 
source, smeared-sink correlators to extract the ground- 
state energy differences and to handle the associated er¬ 
ror correlations. One such example is the spin-averaged 
IP — IS” mass difference: 


AMps = [5E{xb2) + 3E{xbi) + E{xbo)]/9 - 

[3E{T) + E{fj,)]/4 . (7) 


It is this quantity which we use to set the scale for the 
lattices. Figure displays the results versus the pion 
mass. The leftmost point is the chirally extrapolated 24^, 
Nf = 2 + 1 difference (for amq = 1.5 « arub, see below) 
and the closest black point is that for the 32^, Nf = 2 en¬ 
semble (a). Using the bb experimental value of 457 MeV 
[ST], leads to a"! = 2726(206) MeV and 2837(55) MeV, 
respectively. Results from twice the heavy-quark mass 
(amg = 3.0) are also displayed, showing the relatively 
small dependence of this splitting on mg. 

Using ground-state energy levels from the finite- 
momentum QQ vector correlators, we fit the dispersion 
relation to the following form: 


E — Eq + 


‘2Mkir, 




( 8 ) 


Resulting values for (in units of the IP — l^ split¬ 

ting, AMps) as a function of the pion mass are pre¬ 
sented in Fig. The larger error bars on the chirally 
extrapolated 2 U result and on the 32^ result are those 
which also include the error in the lattice spacing de¬ 
termination (about 7.5% and 2%, respectively; the same 
applies to all following figures). The experimental value 
of Mp / AM PS is also plotted and one can see agreement 
when amg = 1.5 on the 24^, Nf = 2 -|- 1 lattices. For 
the 32^, Nf = 2 lattice, a slightly lower value for the 
heavy-quark mass {amg « 1.3) may have been appro¬ 
priate. With the lack of a chiral extrapolation and the 
quenching of the strange quark, however, it is difficult 
to determine which systematics would become absorbed 
into such an adjustment. We use amg = 1.5 as our 
“working value” of the physical bottom-quark mass on 
all ensembles. 

In order to find the radiative corrections for the spin- 
dependent terms (ca and C 4 ) in the heavy-quark action, 
we look at the “spin-orbit” and “tensor” energies of the 
IP bb system: 


Eso = [-2P(x&o) - 3P(x&i) -b 5 P(x&2 )] /9 , (9) 


Et = [-2P(xbo) + 3P(xbi) - E{xb 2 )] /9 . (10) 

These are roughly proportional to C 3 and c|, respectively 
(the corresponding experimental values are 18.20 MeV 



(M^AMpg)® 


FIG. 2: Ground-state vector QQ kinetic mass versus pion 
mass. 


and 5.25 MeV). In Fig. we plot our values for these 
energies for the cases where C 3 = C 4 = 1 and C 3 = 1 , 
C 4 = 1.2. Within the errors the tree-level C 3 = 1 appears 
to work fine for Esoy whereas the correction C 4 = 1.2 
leads to better agreement for Ep- For the Nf = 2 re¬ 
sults, a slightly higher C 4 appears to be needed (at least 
without a chiral extrapolation); together with a lower 
value for heavy-quark mass (amg = 1.3 for better M^m 
and Eso as well), the value C 4 = 1.23 would bring the 
tensor energy in better agreement with experiment. Now 
that the lattice scale and parameters of the heavy-quark 
action have been handled, we can turn our attention to 
other results. 

Another quantity which follows from a jackknife ratio 
analysis is the spin splitting in the ground-state bb S- 
waves: M{T) — M{rj^). In Fig.one can see that the 24^, 
Nf = 2 - 1-1 results extrapolate to a value slightly below 
the experimental value. However, when we include the 
error in the lattice spacing, we see that this discrepancy is 
only a little more than Icr. The Nf = 2 result is also low, 
even after naively extrapolating to amg = 1.3, C4 = 1.23 
(using the amg = 3.0 and C4 = 1 results). 

We also use jackknife ratios to look more closely at 
the IP spin splittings, including the hf,, and the energy 
differences to (and among) the \D states. These results 
for the 24^, Nf = 2-1-1 chiral extrapolation, the Nf = 
2 - 1-1 ensemble closest to the chiral limit (b), and the 32^, 
Nf = 2 ensemble (a) can be found in Table 

As eluded to earlier, we also fit the heavy-quarkonia 
correlator matrices to the form of Eq. H- Depending 
upon the quantum numbers and time interval being con¬ 
sidered, we fit anywhere from one to three energy levels. 


VI 










5 


0.06 




< 

i- 

W 


0.04 


[] 


0.02 


0.00 


0.0 


+ 

w 


(E/AMp^)j,^p 

amn=1.5, c. 


= 1 . 


24'" 

3 


amq=1.5, c^-l, 24 
amq=1.5, c^=l.2, 32'^ 
arrin^l.S, c.^l, 32® 


© 

m 


0.5 


1.0 


(MyAMpg) 


FIG. 3: Spin-orbit (SO) and tensor (T) energies from IP 
QQ states versus pion mass, showing that ca = 1, C 4 = 1.2 
works best. 
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FIG. 4: T — 771 , mass differences versus pion mass. 


This provides us with estimates of the masses of first- 
excited states. The Nf = 2 and the chirally extrapolated 
Nf = 2 + 1 results are compiled in Table |VI| One must 
be careful in interpreting the results of such states. The 
continuum-spin identification in Table [lIT| is the minimum 
possible (Jmin) for the associated irreducible representa¬ 
tion. Ideally, one should create the states of interest with 
operators in different lattice representations and look at 
the amplitudes as a way to confirm the associated spin 


(see, e.g., [351 [SH])- For the Ai operators, after 0, the 
next lowest possible continuum spin is 4. Such states 
should be much higher in mass and we believe that we 
can say with confidence that the first-excited states we 
see for such operators are still spin 0. For the other rep¬ 
resentations, the separation in possible J values is not so 
large: the next lowest values are 3 for Ti and T 2 and 4 
for E. Luckily, for the case of the T 2 operator we have 
the E irreducible representation as well and can verify 
that the first excitation is consistent for both cases (in 
the end we average them). Otherwise, since the possible 
continuum-spin separation for the other operators is at 
least 2 and we are only dealing with first-excited states, 
we assume that J = Jmin for these states as well. 


B. B, Bs, Be mesons 

Just as in the bottomonia case, we create jackknifed 
ratios of smeared-source, smeared-sink correlators to look 
at ground-state mass splittings of B mesons. Table |VII| 
shows the results for the Nf =2-1-1 chiral limit, the 
Nf = 2 + 1 ensemble closest to the chiral limit (b), and 
the Nf = 2 ensemble (a). 

The spin splittings for the S- and P-wave B mesons are 
shown in Figs. and [^ as a function of the pion mass. 

On all ensembles, the B* — B difference agrees with 
experiment. Replacing the light quark with a strange or 
charm one leads to our results for B^ and Be mesons. 
Nf = 2-1-1 ground-state splittings, as well as differences 
to some first-excited states are shown in Table IVTITl One 
may note the fact that the B* — B splitting appears to 
be slightly larger than that for B* — Bg. This may be a 
sign of the small volumes’ effect on the light quarks. 


TABLE VI: Results for bottomonia mass splittings (in MeV) 
using arriQ = 1.5 and C 4 = 1.2 for the Nf = 2-1-1 chiral limit 
(x)i Nf = 2 + 1 ensemble b, and the Nf = 2 ensemble a. 
The first error comes from the fit, the second from the scale 
setting (AMps). 


splitting 

fR — XbO 
iR - Xbi 
Xb2 - iR 
IR - hb 
Vb2 ~ I'S' 
Tm - IS' 

Ab2 - rib2 

n'b- Vb 

x'bo ~ XbO 
Xbl ~ Xbl 
Xb2 ~ Xb2 
K - hb 
Vb2 ~ Vb2 
'^b2 — A b2 


V/ = 2 + 1 (x) V/ = 2 + 1 (b) 


59.7(2.8)(4.5) 
46.8(7.5)(3.5) 
10.4(4.4)(0.8) 
15.4(4.0)(1.2) 
1.3(1.8)(0.1) 
670(150)(50) 
720(190)(55) 
130(95)(10) 
676(95)(51) 
671(97) (51) 
560(190)(40) 
690(200)(50) 
510(220)(40) 
620(200)(50) 
590(350)(45) 
1550(400)(120) 


58.7(1.1)(1.8) 
43.9(3.1)(1.3) 
9.4(1.8)(0.3) 
14.4(1.6)(0.4) 
2.13(64)(6) 
751(58)(23) 
718(77)(22) 
60(40)(5) 
596(31)(18) 
584(33)(18) 
625(69)(19) 
684(76)(21) 
583(86)(18) 
630(79)(19) 
810(170)(30) 
1170(140)(40) 


Nf = 2 (a) 
50.76(90) (98) 
32.9(1.5)(0.6) 
9.04(86)(18) 
12.01(78)(23) 
1.09(50)(2) 
788(35)(15) 
824(47)(16) 
21(18)(1) 
485(36) (9) 
474(38) (9) 
668(93)(13) 
729(97)(14) 
647(78)(13) 
654(79)(13) 
903(66)(18) 
990(190)(20) 
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FIG. 5: B* — B mass differences versus pion mass. 


FIG. 6: B* — Bq mass differences versus pion mass. 


Figure [^displays the B*q — B mass difference as a func- 


TABLE VII: Results for B-meson mass splittings (in MeV) 
using arriQ = 1.5 and C 4 = 1.2 for the Nf =2 + 1 chiral limit 
(x)) =2 + 1 ensemble b, and the Nf = 2 ensemble a. 

The first error comes from the fit, the second from the scale 
setting (AMps). 

splitting Nf =2 + 1 (y) Nf =2 + 1 (b) Nf = 2 (a) 

B* - B 45(12)(3) 43.9(4.5)(1.3) 49.4(3.0)(1.0) 

Bl-Bo 54(11)(4) 46.2(4.1)(1.4) 47.3(3.2)(0.9) 

Bo-B 281(47)(21) 265(25)(8) 357(29)(7) 

B*-B* 253(46)(19) 256(25)(8) 357(29)(7) 


TABLE VIII: Results for B(s_c)-meson mass splittings (in 
MeV) using amg = 1.5 and C 4 = 1.2. for the Nf =2 + 1 
chiral limit (y) and Nf = 2 + 1 ensemble b. The first error 
comes from the fit, the second from the scale setting (AMps). 


splitting 

b; - Bs 

— ^sO 

Bt - B* 
Bto - B 
B*i - B* 
B* - Be 

-^cl ~ -DcO 

B* -B* 
B^o - B 
B,*! - B* 
B'- B, 

TD*' JD* 

-^c -^c 


iV/ = 2 + 1 (x) V/ = 2 + 1 (b) 
43(10)(3) 42.7(3.8)(1.3) 

49.4(9.4)(3.7) 45.4(3.4)(1.4) 

81(21)(6) 30.9(3.1)(0.9) 

418(29)(32) 345(12)(10) 

396(29)(30) 332(13)(10) 

52.0(4.5)(3.9) 55.4(1.4)(1.7) 

46(20)(3) 56.5(6.5)(1.7) 

1165(79)(88) 1119(28)(34) 

1405(79)(106) 1417(31)(43) 

1346(82)(102) 1404(31)(43) 

560(260)(40) 780(120)(25) 

600(250)(45) 780(100)(25) 


Bio-Bio 460(360)(35) 750(120)(25) 


Bl^-Bl^ 460(310)(35) 796(120)(25) 


0.2 


m 


PQ 


0.1 


0.0 


0.00 


O amq=1.5 + aM 

O amr,=3.0 


K.exp 


0.01 


0.02 
(a MJ = 


0.03 


0.04 


FIG. 7: Bio — B mass differences versus pion mass. Indica¬ 
tions are that the B*o is below the BK threshold. 


tion of pion (and thereby the kaon) mass. All indications 
are that the i?*Q is below the BK threshold. The same is 
true for the B^ and the B*K threshold (see Table VIIII. 

Due to the rather large amount of smearing used for 
the light- and strange-quark sources on ensembles b and 
e [S3] (see also Table 0, we are not able to resolve radi¬ 


ally excited B^*'> or Bs' mesons on these ensembles (and 
therefore, not in the chiral limit either). For the present 
study, light- and strange-quark sources with much less 
smearing were created on ensembles c and d to better 
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m 


< 


1.5 


1.0 — 


0.5 — 


0.0 


0.0 


® (AM/AMpg),^p 

C4 = 

w amq=1.5, € 4 =!.2, 32 


0.5 


1.0 


(M^AMpg) 


FIG. 8 : Pseudoscalar B' — Be mass differences versus pion 

mass (775 — rjb and — Bs energies are also plotted for com- FIG. 9: Kb — B mass differences versus M^. 

parison). 


study such states. The results for the ensemble closer to 
the chiral limit (c) are shown in Table The results for 


the corresponding pseudoscalar splitting appear 

in Fig. along with the and 77 ^ — iji, differences. 

Clearly, many more statistics are needed here, especially 
to reach a more reliable chiral limit. 


C. Singly heavy baryons 

The naming conventions for heavy baryons may not 
be strictly obeyed here, especially where negative-parity 


TABLE IX: Results for B(s^c)-meson first-excited-ground 
state mass splittings (in MeV) using amq = 1.5 and C 4 = 1.2 
on ensemble c. The first error comes from the fit, the second 
from the scale setting (AMps). 


splitting 

Nf = 2 + 1 {c) 

B' - B 

670(120)(30) 

B*' - B* 

640(90)(25) 

B*o - Bo* 

790(160)(30) 

Bi*' - Bi* 

760(110)(30) 

B'g-Bg 

697(100)(27) 

Bg — B* 

659(80)(25) 

Bto ~ BgQ 

730(110)(30) 

Bti — B*i 

730(95)(30) 

Be- Be 

661(73)(25) 

B*e - B*e 

662(73)(25) 

n*' D* 

-L>c0 ~ -^cO 

630(110)(25) 

Bh — Bel 

675(95)(26) 


states are concerned. Therefore, for the purpose of nota- 
tional clarity, we briefly point out the quantum numbers 
and the corresponding names of the baryon states dis¬ 
cussed in this subsection: 


fife 


= i + (0) 


= i + (l) 

^b , ^ 

= f + (l) 

V* 77* 
^b ’ ^b 

J^{s,,) = r{0) 

A* 

1 ^b 

J^{s,,) = r{l) 

^b ’ “6 




where 


^dq 


IS 


the spin of the diquark appearing in the 
associated baryon interpolator {not that of the state). 

Table 1^ lists the results for our mass splittings among 
the singly heavy baryons. 

Figurej^shows the A^ — B mass difference as a function 
of The results on all lattices are high when compared 
to experiment, a sign that the small volumes (see Table 
[I]) may be drastically affecting our results for baryons 
with two light quarks. The Kb — B splitting appears 
again in Fig. along with the the analogous — Bg 
splitting, but here again our lattice results are too high, 
536(78) MeV (y-extrap. Nf = 2-1-1), when compared 
to experiment (427 MeV). The fib — Bg difference shows 
much better agreement. 

Mass differences from alternative spin-flavor combi¬ 
nations appear in Fig. 11 Here, the (spin averaged) 


S), — Ab and EJ — Sb splittings also show possible signs 
of finite-volume-induced enhancements (the dotted re¬ 
sults are from chiral extrapolations without the heaviest 
pion mass, e). Other spin and multiplet splittings can be 
found in Table (sometimes also with the same limited 
chiral extrapolation). 

Figure displays the mass differences between the 
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FIG. 12: Mass differences - involving the and Ei,i, 

FIG. 10: =.bb — T, ^bc — Be, ^b — Bs, and Kb — B mass negative-parity states - versus M^. 
differences versus M^. 


D. Doubly heavy baryons 


< 


w 


0.8 


0.6 — 


0.4 


w 0.2 


0.0 


0.0 


^ (AM/AMpg),,p 

□ amQ=1.5, 04 = 1 . 2 , 24'" 
K amQ=1.5, 04 = 1 . 2 , 32" 


0 i 


0.5 1.0 


(M^AMpg) 


FIG. 11: — Kb and EJ — Et mass differences versus M^. 


Kb and its parity partner, A^. Again, there is a clear 
overestimate here; the experimental value being around 
300 MeV [20]. Mass differences involving this and other 
negative-parity states can be found in Table Em- 


As before, we point out the names of the states dis¬ 
cussed in this subsection: 

^bb 


= i + (0) : 


= ^(1) : 


■^bc 


7' O' " 
^bc ^ ^‘‘bc ’ ^bb 


TABLE X: Results for bqq' (where q, q' = u, d, a) baryon 
mass splittings (in MeV) using amg = 1.5 and a = 1.2 for 
the N[ =2-1-1 chiral limit (y), the Nf = 2 -|- 1 ensemble b, 
and the Nf = 2 ensemble a. The first error comes from the 
fit, the second from the scale setting (AMps). (^Using only 
the three lightest values.) 


splitting 

A/ = 2 + 1 (y) A/ = 2 + 1 (b) 

Nf = 2 (a) 

Kb-B 

Eb — Bs 
Q,b — Bs 

450(99)(34) 
536(78) (40) 
746(43)(56) 

474(26)(14) 

528(20)(16) 

695(16)(21) 

444(12)(9) 

Hi — Kb 

231(26)(17) 

87.2(8.2)(2.6) 

- 

E‘‘^ - Ai 

322(62)(24)t 

274(19)(8) 

229(14) (4) 

“6 “ “fa 

238(25)(18) 

190(10)(6) 

- 


154(40)(12)t 




314(34)(24) 

221(12)(7) 

- 

Si* - Ei 

32(19)(2)f 

31.9(5.1)(1.0) 

31.9(3.0)(0.6) 

1 — 1 * 1 — 1 / 

36.6(9.8)(2.8) 

33.4(3.1)(1.0) 

- 

fii - fii 

31.6(8.5)(2.4) 

27.6(3.2)(0.8) 

- 

Ai* - Ai 

514(67) (39) 

455(26)(14) 

495(27)(10) 

' 

—b — 

385(79)(29)t 

398(20)(12) 

- 

Ei* - Ei 

142(30)(11) 

183(12)(6) 

274(15)(5) 

Sr - Ai* 

53(14)(4) 

46.1(6.4)(1.4) 

36.0(4.9)(0.7) 

1 — 1 /* t— 1 :4; 

“b ~ “b 

262(36)(20) 

237(15)(7) 

- 

n': - Qt 

308(28)(23) 

249(11)(8) 

- 
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= 1^(1) 

^^(s.,) = r(o) 


j^(s,,) = f-(i) 


77* 

'~‘bc 

77" 

^bc 

^bb 

77 '* 

'‘be 


o* 

“be ’ 

Kc 


77 * 

^bb 


77 '* 

-66 




bb 


O'* 

“hb 


and any “radial” excitations are denoted by a (2) after 
the name. 

Mass differences involving the 6c-baryons can be found 
in Table pCT] Those for doubly bottom baryons appear in 
Table IXlil 

Also present in the previously mentioned Fig. are 
the results for the (yet to be observed) Etc — Be mass 
difference. Using the physical value for M{Bc), we find: 
M(Ebc) = 6887(103)(46) MeV (y-extrap. Nf=2 + 1). 

The Sbb — T and — T energies may also be seen in 
Fig. With the physical value for M(T) as input, we 
find: 

M{Ebb) = 10201(10)(14) MeV {Nf = 2) ; 

MiEbb) = 10267(44) (61) MeV (x-extrap. Nf = 2 + 1) ] 
Miribb) = 10356(34)(68) MeV (y-extrap. Nf = 2 + 1). 

Along with the A^, the negative-parity also appears 
in Fig. 

For more results involving spin splittings, negative- 
parity states, and (for the 66-baryons) first-excited states, 
see Tables XI and uniii]. 


E. Triply heavy baryons 

As before, we point out the names of the states dis¬ 
cussed in this subsection: 

J^{Sdq) = 5'''(1) : Ubec , Ubbe 

J^{sdq) = f (1) : ^bcc 1 ^bbe > ^bbb 

J^{sdq) = 5 ( 1 ) : ^'bbe 

J^{sdq) = f (1) : ^bcc j ^'bbe 
and any “radial” excitations are denoted by a (2) after 
the name. 

Mass differences involving bcc and bbc baryons can be 


TABLE XL Results for beq (where q = u,d, s) baryon mass 
splittings (in MeV) using anig = 1.5 and a = 1.2 for the 
Nf = 2-1-1 chiral limit (x) and the Nf = 2-1-1 ensemble b. 
The first error comes from the ht, the second from the scale 
setting (AMps). 


splitting 

Nf = 2 + 1 ix) Nf = 2 + 1 {h) 

“be Be 

611(103)(46) 

681(41)(21) 

^bc “be 

152(34)(11) 

51(13)(2) 

“be “ “fac 
^bc ^bc 

50(19)(4) 
38(15)(3) 

47.2(7.2)(1.4) 

43.5(4.9)(1.3) 

“be '—^bc 

flic “ flic 

26(13)(2) 

21 ( 11 )( 2 ) 

26.7(4.6)(0.8) 

25.0(3.4)(0.8) 

^bc “6c 

^bc “ ^bc 

290(66)(22) 
342(53)(26) 

314(21)(10) 

334(14)(10) 

“be '— 'be 

O'* 0 " 

^^bc “be 

57(17) (4) 
52(14)(4) 

48.1(5.4)(1.5) 

45.6(4.0)(1.4) 


0.10 


<i 0.05 


0.00 


0.00 



m(j=1.5, C4=1.2. 24“ 

□ M(n,J-2M(Bj + M(T)/2 
o M(n,,J-M(Bj-M(T)/2 
X M(n,,J-3M(T)/2 


0.01 


0.02 

(aMj^ 


0.03 


0.04 


FIG. 13: Slice, Slbbc, and Qbbb masses versus M^. 


found in Table XIII Triply bottom results appear in Ta¬ 
ble EVl 

With the physical values for M{Be) and M(T) as in¬ 
put, we find: 

M{nbcc) = 7984(27)(12) MeV (x-extrap. Nf = 2 + 1) ; 
M(Ubbc) = 11182(27)(13) MeV (x-extrap. Nf = 2 + 1). 


TABLE XII: Results for bbq (where q = u,d, s) baryon mass 
splittings (in MeV) using amg = 1.5 and C 4 = 1.2 for the 
Nf = 2-1-1 chiral limit (x), the Nf = 2 -|- 1 ensemble b, and 
the Nf = 2 ensemble a. The first error comes from the fit, 
the second from the scale setting (AMps). (^Using only the 
three lightest M,r values.) 


splitting 

Nf = 2 + 1 [x) Nf = 2 + 1 ih) 

Nf = 2 (a) 

Ebb-r 

807(44)(61) 

802(17)(24) 

741(10)(14) 

Qbb — T 

896(34)(68) 

841(12)(25) 

- 

“bb “66 

49(12)(4) 

34(20)(3)^ 

33.7(4.9)(1.0) 

35.1(6.5)(0.7) 

Qlb — ^bb 

45.1(9.5)(3.4) 

31(17)(2)1' 

32.5(3.7)(1.0) 

■ 


295(41)(22) 

293(17)(9) 

341(17)(7) 

— Qbb 

362(44)(27) 

299(19)(9) 

- 

“bb “ “bb 

54(12)(4) 

40(20)(3)1^ 

40.5(6.9)(1.2) 

28.4(3.4)(0.6) 

^'bb — ^bb 

50(10)(4) 

37(17)(3)1' 

37.3(5.0)(1.1) 

■ 

“bb (2) “bb 

- 

671(69)(20) 

375(110)(7) 

“bb (2) “bb 

- 

591(63)(18) 

387(83) (7) 

^bb(2) ^bb 

- 

624(54)(19) 

- 

n*bb{2)-n*bb 

- 

564(50)(17) 

- 

“66(2) — E'bb 

- 

739(92)(22) 

- 

“66(2) — Abb 

- 

793(95)(24) 

- 

1166(2) — ^'bb 

- 

722(74)(22) 

- 

116*6(2) - ll6*6 

- 

711(71)(22) 

- 
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m 

a. 



31.6 — 


31.4 


31.2 


[| 
— ^:> 


31.0 


0.0 




amq=1.5, 04 =!.2, 24 
n M(nbbb)-3M(T)/2 
o M(T) + M(B*)-M(nbbb) 
amo=1.5, € 4 = 1 . 2 , 32^ 


0.5 1.0 


(M^AMpg)^ 


FIG. 14: Q.bbb masses versus M^. 


With the physical value for M(T) as input, we find: 
MiVLbbb) = 14357.6(4.4)(3.3) MeV {Nj = 2) ; 

M{Vtbbb) = 14369(21)(14) MeV (y-extrap. TVy = 2 + 1). 

Figure 13 shows the flbcc, and 0.bbb masses (ref¬ 
erenced from an appropriate combination of M(T) and 
M{Bc)) versus the pion mass. Once the heavy-quark 
masses are “subtracted”, seemingly little difference re¬ 
mains between these triply heavy systems. 

Figure displays the mass achieved via two 
schemes: from E{V,i,bb) — 3E{T)/2 and M(T); or from 
E{T) + E{B*) - E{nbbb) and M(T) -k M{B*). The two 
estimates agree in the Nf = 2-1-1 chiral limit and on the 
Nf = 2 lattice. 

For further results (e.g., spin and parity splittings) in¬ 
volving hcc and bbc baryons, see Table |XIII| Estimates 
for the first-excited Vlbbb appear in Table |XfV 


Using (jackknifed) combinations of four correlators, we 
were able to compare the effect of trading one of the b 
quarks for a lighter one in the mesons with the same 


effect in the baryons. These results appear in Table [XV] 
For the first meson difference, T — Be, the spin changes 
from 1 to 0, but these states are much better determined 
(experimentally) than the corresponding rjb and B* (this 
extra spin jump adds about 50 MeV to the first four 
rows). Considering the errors, there appears to be little 
difference between the splittings involving T — B* and 
T — B *, but exchanging the b quark for a strange or light 
one appears to have a larger effect on the baryons than 
exchanging 6 —)■ c, especially when looking at baryons 
with fewer heavy quarks [T]. 


IV. DISCUSSION 

We reserve this section for a brief discussion of possible 
sources of systematic errors affecting our results [T]. 

(1) Finite-volume effects: 

As can be seen in Table |Ij all volumes are rather small 
when compared to the dynamical pion masses {Mt^L < 
4), with L ranging from 1.7-1.8 fm for the Nf = 2 + 1 
ensembles to 2.2 fm for the Nf = 2 one. Obvious energy 
enhancements can be seen in some hadrons containing 


TABLE XIV: Results for bbb baryon mass splittings (in MeV) 
using aniQ = 1.5 and C 4 = 1.2 for the Nf = 2-1-1 chiral limit 
(x), the Nf = 2 -|- 1 ensemble b, and the Nf = 2 ensemble a. 
The first error comes from the fit, the second from the scale 
setting (AMps). (^Constant fit, using only the three lightest 
Mtt values.) 


splitting 

Nf =2 + 1 ix) Nf = 2 + 1 {h) 

Nf = 2 (a) 

Hbbb - 3T/2 
T + B*... 

179(21)(14) 

172(10)(5) 

167.6(4.4)(3.3) 

^bbb 

162(31)(12) 

150(13)(5) 

169(12)(3) 

^bbb 

211(128)(16) 
438(33) (33)* 

411(45)(12) 

480(64)(9) 


TABLE XV: Comparison of meson and baryon mass split¬ 
tings (in MeV) resulting from the replacement of one b quark 
with a lighter one. Using arriQ = 1.5 and C 4 = 1.2 for the 
Nf = 2-1-1 chiral limit (x) and the Nf = 2-1-1 ensemble b. 
The first error comes from the fit, the second from the scale 


TABLE XIII: Results for bbc and bcc baryon mass splittings 
(in MeV) using arrig = 1.5 and C 4 = 1.2 for the Nf = 2 + 1 
chiral limit (x) and the Nf = 2 + 1 ensemble b. The first 
error comes from the fit, the second from the scale setting 
(AMps). 


splitting 

A/ = 2 + 1 (x) V/ = 2 + 1 (b) 

^bcc ~ ‘^Bc + T/2 
Q-bbc — Be — T/2 

164(27)(12) 

177(27)(13) 

185(9)(6) 

179(9)(5) 

^bcc ^bcc 

^bbe ~ ^bbe 

28.8(5.6)(2.2) 

21.8(7.0)(1.7) 

27.2(1.8)(0.8) 

26.5(2.2)(0.8) 

O'* 0 * 

^^bcc ^^bcc 

^bbe ~ ^bbe 

449(41) (34) 
375(74)(28) 

422(15)(13) 

372(26)(11) 

TV* rv 

^^bbe ^^bbe 

30(12)(2) 

39.6(4.1)(1.2) 


setting (AMps). 


splitting 


Nf = 2 + 1 (x) Nf = 2 +1(h) 

(T - Be) 

— (Qbbb — l^bbc) 

45(14)(3) 

42.3(3.5)(1.3) 

(T - Be) 

{^bbe ^6cc) 

16.7(9.8)(1.3) 

16.4(3.4)(0.5) 

(T - Be) 

— (Hbb — li'be) 

17(32)(1) 

18(12)(1) 

(T - Be) 

~ i^bb ~ 

34(36) (3) 

27(13)(1) 

(T-Bt) 

— (l^bbb — l^bb) 

-29(20)(2) 

-29.3(8.4)(0.9) 

{T-Bt) 

— (libbc — U(^) 

-91(27)(7) 

-84.8(8.8)(2.6) 

(T-B*) 

^b) 

-130(50)(10) 

-137(17)(4) 

(T-Bt) 

- (Sbb - H() 

-171(65)(13) 

-154(27)(5) 

(T-B*) 

— (libbb — S.bb) 

-35(20)(3) 

-29.7(7.9)(0.9) 

(T-B*) 

~ {^bbe ~ “be) 

-112(30)(8) 

-96(10)(3) 

(T-B*) 

— {Q.bb — -b) 

-137(57)(10) 

-139(21)(4) 

(T-B*) 

— i—bb — Si,) 

-179(74)(13) 

-153(33)(5) 
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light quarks: e.g., A[*\ possibly even 

_B* — B. We can only imagine how large such effects may 
be for higher excitations. However, there are no indi¬ 
cations that hadrons containing strange quarks as their 
lightest constituents suffer from these small boxes. 

(2) Different previous lattice scales; too heavy s and c 
quarks: 

During the tuning and running stages of the present 
ensembles [501 IS], different scales were used than the 
one in the current study. The same was true during a 
tuning of the charm quark [SU [55]. The fact that we 
now use AMps from the bottomonia system to set the 
scale leads to a seemingly ~ 10% enhancement of the c 
and s (in the chiral limit) masses. 

(3) For Nf = 2, too heavy b quark: 

Looking at Fig. one can see that on the Nf = 2 
ensemble, arug = 1.5 is too large (by ^ 10%) for the 
bottom quark. 

(4) Different smearings; different excited-state “con¬ 
tamination” : 

On ensembles c and d, much less smearing was chosen 
for the light and strange quarks in an attempt to better 
identify excited states. However, this may lead to dif¬ 
ferent levels of contamination from excited states when 
considering the ratios of smeared-smeared correlators on 
the different ensembles leading to the Nf = 2 + 1 chiral 
limit. 

(5) Relatively heavy u, d quarks; long chiral extrapo¬ 
lations: 

For many systems considered herein, this one is not 
so much a systematic error, outside of the fact that it 
systematically causes larger errors in the chiral limit. 
More trouble arises from this when the system consid¬ 
ered should cross a threshold between ensemble b and 
the chiral limit. 


V. CONCLUSIONS AND OUTLOOK 

We have presented the heavy-hadron spectrum aris¬ 
ing from NRQCD-approximated b quarks and improved 
clover-Wilson c,s,d,u quarks on Nf = 2 and 2 + 1 lat¬ 
tices. Singly, doubly, and triply heavy (6, c) systems were 
considered and results were found for spin splittings, al¬ 


ternate parities, and in some cases, “radial” excitations. 
Relative mass differences between mesons and baryons 
resulting from the exchange of one b quark for a lighter 
one were presented as well. A number of systematics 
were pointed out (e.g., small volumes; see above), but 
these have not been fully quantified and the reader must 
use some caution with the results. Perhaps more trust¬ 
worthy are the lower-lying spectra (no radial excitations) 
of hadrons containing only 6, c, or s valence quarks. A 
more careful analysis of all the data generated (e.g., with 
anig = 3.0 or the 40^, Nf = 2 correlators) and further 
runs on existing 32^, Nf = 2 + 1 ensembles could uncover 
more about the errors incurred, and an expansion of the 
code [HI] could lead to a study of possible bbqq states, but 
we must excuse ourselves from taking this project further 
as we have precious little free time and have found new 
ways in which to be wrong [69] . 
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